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SPT-CL J0546-5345: A MASSIVE Z > 1 GALAXY CLUSTER SELECTED VIA THE SUNYAEV-ZEL'DOVICH 

EFFECT WITH THE SOUTH POLE TELESCOPE 

M. Brodwin/-^-^^ J. RuEL,^ P. A. R. Ade,'' K. A. Aird,'' K. Andersson," M. L. N.Ashby/ M. Bautz,^ G. Bazin/'* 
B. A. Benson,"-"' L. E. Bleem,^ " J. E. Carlstrom,^'1°'"'1^ C. L. Chang,'-''"' T. M. Crawford,^-!^ A. T. Crites,^-^^ 
T. DE Haan,"* S. Desai," M. a. Dobbs,"* J. P. DuDLEY.i^* G. G. Fazio, ^ R. J. Foley,"^'' W. R. Forman,^ 
G. Garmire,^^ E. M. George,^"* M. D. Gladders,"'^^ A. H. Gonzalez," N. W. Halverson,^*^ F. W. High,^ 
G. P. Holder," W. L. Holzapfel," J. D. Hrubes,-^ C. Jones, ^ M. Joy,i'' R. Keisler,"-" L. Knox,^" A. T. Lee,^''''^^ 

E. M. Leitch,^'^^ M. Lueker,^'' D. P. Marrone,''-" J. J. McMahon,^'!"'^^ J. Mehl,^-!^ S. S. Meyer,'-'-^'''"'" 
J. J. MOHR,^'**'^^'^^ T. E. Montroy,^'' S. S. Murray,^ S. Padin,'^-^^ T. Plagge,"-^" C. Pryke," "'-^^ C. L. Reichardt," 
A. Rest,3 J. E. RuHL,^'* K. K. Schaffer,^ "' L. Shaw,"''-^'^ E. Shirokoff,!" J. Song," H. G. Spieler,^! B. Stalder,^ 
S. A. Stanford,^" Z. Staniszewski,^^ A. A. Stark, ^ C. W. Stubbs,"'^ K. Vanderlinde," J. D. Vieira,'^ " A. 
ViKHLiNiN,^ R. Williamson, "'^^ Y. Yang," O. Zahn,"' and A. Zenteno""'* 

ApJ, in press 

ABSTRACT 

We report the spectroscopic confirmation of SPT-CL J0546-5345 at (z) = 1.067. To date this is 
the most distant cluster to be spectroscopically confirmed from the 2008 South Pole Telescope (SPT) 
catalog, and indeed the first z > 1 cluster discovered by the Sunyaev-Zel'dovich Effect (SZE). We 
identify 21 secure spectroscopic members within 0.9 Mpc of the SPT cluster position, 18 of which 
are quiescent, early-type galaxies. From these quiescent galaxies we obtain a velocity dispersion of 
1179 

-167 km/s, ranking SPT-CL J0546-5345 as the most dynamically massive cluster yet discovered 
at z > 1. Assuming that SPT-CL J0546-5345 is virialized, this implies a dynamical mass of M200 = 
l-Oito'4 X 10^^ Mq, in agreement with the X-ray and SZE mass measurements. Combining masses 
from several independent measures leads to a best-estimate mass of M200 = (7.95 ± 0.92) x 10^'' AIq. 
The spectroscopic confirmation of SPT-CL J0546-5345, discovered in the wide-angle, mass-selected 
SPT cluster survey, marks the onset of the high redshift SZE-selected galaxy cluster era. 
Subject headings: galaxies: clusters: individual (SPT-CL J0546-5345) — galaxies: distances and 
redshifts — galaxies: evolution 
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1. INTRODUCTION 

Formed from the extreme peaks of the primordial den- 
sity field, galaxy clusters are both a sensitive probe of 
cosmology and an excellent laboratory for galaxy evolu- 
tion studies. Although massive galaxy clusters are easy 
to identify in the local Universe, finding such objects at 
large lookback times has been extremely challenging. A 
patchwork of optical, infrared and X-ray techniques has 
led to a heterogeneous sample of distant galaxy clusters 
(e.g., Stanford et al. 1997; Rosati et al. 2004; Gladders 
and Yee 2005; Mullis et al. 2005; Brodwin et al. 2006; 
Elston et al. 2006; Stanford et al. 2006; Eisenhardt et al. 
2008; Muzzin et al. 2009; Papovich et al. 2010). While 
each technique offers unique advantages, namely, obser- 
vational efficiency (optical) , redshift reach and mass sen- 
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sitivity (IR), and a direct mass-observable (X-ray), none 
produce purely mass-selected, nearly redshift indepen- 
dent cluster samples. 

A long-awaited advance is the arrival of Sunyaev- 
Zel'dovich Effect (SZE) cluster surveys. The SZE is a 
small distortion in the Cosmic Microwave Background 
(CMB) spectrum caused by the inverse Compton scat- 
tering of CMB photons with the hot intracluster medium 
(ICM; Sunyaev and Zcl'dovich 1972). The surface bright- 
ness of the SZE does not depend on a cluster's distance, 
which makes it a nearly redshift- independent method for 
finding massive galaxy clusters (e.g., Carlstrom ct al. 
2002). A new generation of dedicated millimeter- wave 
(mm-wave) experiments to search for galaxy clusters are 
currently surveying large areas of sky, including the Ata- 
cama Cosmology Telescope (Fowler et al. 2007) and the 
South Pole Telescope (Carlstrom et al. 2009). 

Taking advantage of the excellent conditions for mm- 
wave observations at the South Pole, the 10-m South 
Pole Telescope (SPT; Carlstrom et al. 2009) is midway 
through a > 2000 deg^ survey sensitive to galaxy clus- 
ters above « 5 X lO^'^ Mq at all rcdshifts. The first 
three clusters discovered by the SZE were reported by 
Staniszewski et al. (2009, hereafter S09). The complete 
5 cr cluster catalog from the 2008 SPT season, a study of 
the optical properties of these clusters, and a description 
of the survey strategy and goals are given in Vanderlinde 
et al. (2010^ hereafter VIO) and High et al. (2010, here- 
after HIO). 

Here we report the spectroscopic confirmation, velocity 
dispersion and dynamical mass of the first z > 1 SZE- 

selected cluster, SPT-CL J0546-5345, at (z) = 1.067. In 
converting from rsoo to r2oo and A'/5oo to Af2oo through- 
out, we assume an NFW density profile and the mass- 
concentration relation of Duffy et al. (2008). Unless oth- 
erwise indicated, we use Vega magnitudes and adopt a 
WMAP5 cosmology (Dunklcy et al. 2009), with VIm = 
0.264, Ha = 0.736, as = 0.80 and Hq = 71 km/s/Mpc. 

2. OBSERVATIONS 

2.1. SPT-SZE 

In 2008, the SPT survey observed a 196 deg^ area (178 
deg^ after point-source and non-uniform coverage mask- 
ing) for ~ 1500 hr (VIO), reaching a 5 cr point source 
sensitivity of ~ 7 mJy at 150 GHz (2.0 mm). Details of 
the data processing, map-making, cluster extraction and 
significance estimates arc given in S09 and VIO. SPT- 
CL J0546-5345, first reported in S09, has an SZE decre- 
ment of S/N = 7.69 at 150 GHz. It is among the most 
significant of the 2008 cluster candidates with a ^ > 1 
optical red-sequence (HIO). SZE significance contours at 
150 GHz are shown in Figure 1 (Ze/t). 

2.2. Spitzer/IRAC and Optical Imaging 

Mid-infrared Spitzer/IRAC imaging was obtained in 
2009 September as part of a larger program to follow 

up clusters identified in the SPT survey (PID 60099; 
PI Brodwin). The on-target observations consisted of 
8 X 100 s and 6x30 s dithered exposures at 3.6 and 4.5 /im, 
respectively. The deep 3.6 /xm observations are sensitive 
to passively-evolving cluster galaxies down to 0.1 L* at 
z = 1.5. The data were reduced following the method of 
Ashby et al. (2009). Briefly, we correct for column pull- 
down, mosaic the individual exposures, resample to 0'.'86 



pixels (half the solid angle of the native IRAC pixels), 
and reject cosmic rays. 

Optical data were provided by the Blanco Cosmology 
Survey (BCS), an NOAO survey program with deep griz 
imaging over 100 deg^ in a subset of the SPT 05hr and 
23hr fields (Ngeow et al. 2009). The Stellar Locus Re- 
gression method of High et al. (2009) was used for photo- 
metric calibration, extinction correction and Galactic de- 
reddening. The independent calibration from the DES 
data management system (Mohr et al. 2008) gave con- 
sistent results. The BCS data typically reach 5 cr AB 
depths of 24.75, 24.65, 24.35, and 23.5 in griz. Sec S09 
and HIO for a more detailed description of the optical 
data. 

Along with the r and i BCS photometry, the IRAC 
3.6 /xm data were used to select promising cluster mem- 
ber candidates for spectroscopy as described below. A 
false-color optical (ri) + IRAC (3.6 fxm) image is shown 
in Figure 1 (right). The optical-only (grz) image (Figure 
1, left) indicates how optically (rest- frame UV) faint the 
cluster members are, and highlights the importance of 
Spitzer/IRAC imaging for high-redshift cluster detection 
and study. 

2.3. Chandra 

SPT-CL J0546-5345 was observed with 
Chandra/ AClS-l on four separate occasions, for a 
total exposm-c time of 55.6 ks (PID 09800046: PI 
Garmire, PID 11800471; PI Mohr). This yielded 
1304 source counts in the 0.5 — 7.0 keV energy 
range within 0.5 rsoo- Background was subtracted 
using blank sky data normalized to the fiux in the 
9.5 — 12 keV energy band. The cluster luminosity is 
Lx(0.5 - 2keV) = '(6.4 ± 0.4) x lO'^^ erg s^^, and the 
temperature, estimated in the (0.15 — 1) rsoo annulus, 
is Tx = 7.5 'tl'l kcV. The gas mass was estimated 
to be Mg,5oo = 7.3l°:3 x lO^'^ Af© from the X-ray 
surface brightness in the 0.5 — 2.0 kcV band following 
the analysis method described in (Vikhlinin et al. 2006), 
where the gas density is assumed to follow a modified 
^-model (Cavaliere and Fusco-Femiano 1978) that is 
allowed to steepen at large radii and to have a power-law 
cusp at the cluster center. The Chandra X-ray contoiirs 
are overlaid in Figure 1 (right). The reader is referred 
to Andersson et al. (2010, hereafter AlO) for a full 
description of the X-ray analysis of SPT-CL J0546-5345, 
as well as a study of the X-ray properties of 15 of the 
SZE-selected clusters from the VIO sample. 

2.4. Optical Spectroscopy 

Multislit spectroscopic observations were acquired on 
the 6.5-meter Magellan Baade telescope on UT 2010 
February 11. We used the Gladders Image-Slicing Multi- 
Slit Option (GISMO^*, Gladders et al. in prep) mod- 
ule on the Inamori Magellan Areal Camera and Spectro- 
graph (IMACS; Dressier et al. 2006). GISMO optically 
remaps the central region of the IMACS field-of-view 
(roughly 3.5' x 3.2') to sixteen evenly-spaced regions of 
the focal plane, allowing for a large density of slitlcts in 
the cluster core while minimizing trace overlaps on the 

http;//www.lco.cl/telescopes-information/magellan/ 
instruments / imacs /gismo 
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Fig. 1.— Left: Optical 4' x 4' color image (grz) of SPT-CL J0546-5345, with SZE significance contours overlaid (S/N = 2, 4, and 6). 
Right: False color optical (■rj)-|-IRAC (3.6 /xm) image of SPT-CL J0546-5345, with Chandra X-ray contours overlaid (0.25, 0.4, 0.85 and 
1.6 counts per 2" X 2" pixel per 55.6 ks in the 0.5-2 keV band). North is up, east is to the left. Due to its high angular resolution, Chandra 
is able to resolve substructure to the SW, which may be evidence of a possible merger. These images highlight the importance of IRAC 
imaging in studying the galaxies in high redshift, optically faint clusters. Spectroscopic early- type (late- type) members are indicated with 
yellow (cyan) circles. Green squares show the spectroscopic non-members. 



CCD. We used the f/2 camera for its greater red sensi- 
tivity, the 300 l/mm "red" grism and the WB6300-9500 
filter. The seeing was excellent (< 0.5") throughout the 
six 30 min exposures. 

In designing the multislit mask, galaxies were assigned 
a weight proportional to their i-band brightness and their 
proximity to the cluster center, and inversely propor- 
tional to their distance in color space, in both r — i and 
i — [3.6], from the predicted color from a Bruzual and 
Chariot (2003, hereafter BC03) passively evolving 100 
Myr burst model with formation redshift = 3. 

The COSMOS reduction package^^ was used for stan- 
dard CCD processing, resulting in wavelength-calibrated 
2D spectra. The ID spectra were then extracted from 
the sum of the reduced data. We flux calibrated the data 
and removed telluric absorption using the continuum of 
a spectrophotometric standard (Foley et al. 2003). 

3. RESULTS 
3.1. Redshifts 

Redshifts for galaxy spectra exhibiting significant ab- 
sorption features were found through cross correlation 
using the RVSAO package in IRAF (Kurtz and Mink 
1998). The validity of the cross-correlation redshift was 
checked by visual inspection and judged by the presence 
of visible absorption lines; most of these displayed clear 
Ca HK lines. For three members the redshift was deter- 
mined from a strong, unresolved [Oil] emission line. Red- 
shift uncertainties were estimated as twice those given by 
RVSAO (Quintana et al. 2000). An independent extrac- 
tion and visual redshift determination yielded identical 
redshifts for all cluster members within the errors. 

From these Magellan observations we obtained 28 se- 
cure redshifts, of which 21 are spectroscopic cluster mem- 
bers within r < 0.9 Mpc of the SPT cluster position. 
This high success rate validates our spectroscopic selec- 

http:/ /obs. carnegiescience.edu/Code/cosmos/ 



tion algorithm. Representative spectra are shown in Fig- 
ure 2. The majority of the members (18) are quiescent, 
with redshifts determined from Ca HK absorption, as 
expected given their red-sequence selection, central loca- 
tion and the extreme mass of the cluster. 

Rest Wavelength (A) 
3200 3400 3600 3800 4000 4200 




'I I 

6500 7000 7500 8000 8500 



Observed Wavelength (A) 

Fig. 2. — Representative spectra of cluster member galaxies for 
SPT-CL J0546-5345. Vertical dotted lines indicate [Oil] and Ca 
HK features. Most (18/21) confirmed members are passive, early- 
type galaxies. 



3.2. Velocity Dispersion 

An iterative 3 a clipping algorithm was used to identify 
cluster members. We list the 21 secure cluster members 
in Table 1, and plot the redshift histogram in Figure 3 
(left). In this Figure the 18 passive members, for which 
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TABLE 1 

Spectroscopic members of SPT-CL J0546-5345 





R.A. 


Dec. 




Principal 


ID 


(J2000) 


(J2000) 


z 


(52^ Spectral Feature 



J054637.2- 
J054636.2- 
J054642.3- 
J054633.6- 
J054635.1- 
J054638.9- 
J054634.8- 
J054635.8- 
J054638.6- 
J054644.3- 
J054645.0- 
J054641.8- 
J054630.3- 
J054638.8- 
J054633.8- 
J054646.8- 
J054634.9- 
J054636.4- 
J054644.4- 
J054646.8- 
J054642.0- 



534525 
534413 
534554 
534540 
534502 
534356 
534555 
534542 
534513 
534538 
534625 
534613 
534425 
534620 
534409 
534610 
534437 
534619 
534502 
534557 
534543 



05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 
05:46: 



37.23 
36.22 
42.36 
33.67 
35.14 
38.90 
34.81 
35.82 
38.68 
44.36 
45.04 
41.84 
30.37 
38.80 
33.86 
46.89 
34.92 
36.43 
44.42 
46.87 
42.08 



-53:45:25.6 
-53:44:13.4 
-53:45:54.5 
-53:45:40.5 
-53:45:02.0 
-53:43:56.2 
-53:45:55.6 
-53:45:42.0 
-53:45:13.7 
-53:45:38.4 
-53:46:25.9 
-53:46:13.5 
-53:44:25.3 
-53:46:20.2 
-53:44:09.3 
-53:46:10.9 
-53:44:37.5 
-53:46:19.5 
-53:45:02.7 
-53:45:57.6 
-53:45:43.9 



1.0647 
1.0656 
1.0548 
1.0775 
1.0676 
1.0743 
1.0567 
1.0710 
1.0692 
1.0619 
1.0662 
1.0712 
1.0502 
1.0681 
1.0586 
1.0625 
1.0705 
1.0805 
1.0661 
1.0686 
1.0524 



0.0002 
0.0004 
0.0003 
0.0005 
0.0003 
0.0003 
0.0002 
0.0003 
0.0003 
0.0002 
0.0003 
0.0004 
0.0005 
0.0004 
0.0005 
0.0005 
0.0005 
0.0007 
0.0025 
0.0002 
0.0002 



Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 
Ca HK/D4000 

on 
on 
on 



^ Redshift errors are twice those given by RVSAO. 




1.02 1.04 1.06 1.08 1.10 1.12 
Redshift 



Fig. 3. — Le/t: Histogram of spectroscopic redshifts obtained for SPT-CL J0546-5345, showing a very broad peak consisting of 21 
members (shaded regions) centered at (z) = 1.0665 (dotted line). The passive galaxies used in the dynamical analysis are shaded red. 
Galaxies for which redshifts were determined from [Oil] emission are shaded blue. Spectroscopic non-members are unshaded. The inset 
shows all the secure redshifts obtained toward this cluster. Right: Color-magnitude diagram in i — [3.6] within a 2.5' radius around 
SPT-CL J0546-5345. The red (blue) squares are the spectroscopically confirmed passive ([Oil]) members. The passive members form a 
rich red-sequence. Spectroscopic non-members are marked with X symbols. The 5 cr color-mag limit is indicated by the dashed line. One 
member is omitted due to low S/N photometry. The red line is the rest-frame Coma CMR, in the approximately [U — H) color probed by 
our filters, normalized to match a BC03 passively evolving L* galaxy. 



the redshift was obtained from Ca HK absorption hnes, 
are shaded red. Redshifts of the remaining three mem- 
bers, shaded blue, were secured from strong [Oil] emis- 
sion. Spectroscopic non-members are unshaded. From 
the spectroscopic redshift histogram, SPT-CL J0546- 
5345 appears to be a single massive halo, with no evi- 
dence of substructure or merger activity. We note, how- 
ever, that the X-ray image (Figure 1) shows a substruc- 
ture extending to the SW indicating that a minor merger 
may be taking place (AID). 



We use the robust biweight estimator of Beers et al. 
(1990) to estimate the mean velocity and velocity dis- 
persion, applying the relativistic correction and the usual 
correction for velocity errors (Danese et al. 1980). For the 
full membership, we find a mean redshift and velocity dis- 
persion of z = 1.0661 lo:oo22 and a = 1181 +111 km/s, re- 
spectively. The errors in both quantities, obtained from 
bootstrap resampling, represent the 68% confidence in- 
terval. 

Since late- type members are often inf ailing, they tend 
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to yield broader dispersions than early- types (e.g., Gi- 
rardi et al. 1996; Fadda ct al. 1996; Mohr et al. 1996; 
Koranyi and Geller 2000). We therefore explore lim- 
iting our analysis to the 18 early-type members, as 
these are expected to better reflect the cluster poten- 
tial. This approach yields nearly identical results, with 
a mean redshift oi z = 1.0665 ^q qq2i ^^'^ dispersion of 
a = 1179 l^gy km/s. Although the biweight estima- 
tor is optimal, we also compute dispersions for these 18 
members using both the gappcr method (Beers ct al. 
1990) and the simple standard deviation. These yield 
fjgap = 1170^128 km/s and ctsd = 1138 1^^^ km/s, re- 
spectively, both in excellent agreement with the biweight 
dispersion. There is no evidence of merger activity in 
SPT-CL J0546-5345 from the line-of-sight velocities, in 
either the appearance of Figure 3 or in the dispersion 
measures, although our redshift sampling is too sparse 
to rule it out. 

3.3. Color-Magnitude Relation 

We plot the i — [3.6] color-magnitude diagram for galax- 
ies within a 2.5' (1.22 Mpc) radius of SPT-CL J0546-5345 
in Figure 3 (right). The spcctroscopically confirmed pas- 
sive galaxies define a rich red sequence, indicative of a 
cluster that is already well-evolved at z = 1.067. The 
dashed line indicates the 5a color-mag limit; one con- 
firmed member is not detected robustly enough in the 
optical for inclusion in this plot. The solid red line indi- 
cates the slope of the rest- frame color- magnitude relation 
(CMR) of Coma (Eisenhardt et al. 2007), normalized to 
L* at 3.6 /xm in the passively evolving Zf = 3 BC03 model 
described above. 

3.4. Dynamical Mass 

There are several approaches in the literature for cal- 
culating the dynamical mass. Perhaps the most widely 
adopted is that of Carlberg et al. (1997), who use the 
simple definition M200 = 4/3 7rr2oo 200pc, where r2oo = 
V3a/10H{z) and Pc is the critical density. Although 
this approach provides a convenient basis for comparison 
with other studies, it relics on the simplifying assumption 
that clusters are singular isothermal spheres. We there- 
fore prefer to use the simulation-based a — M200 relation 
of Evrard et al. (2008), 

h{z) \crDM,15 / 

where (Tdm,15 ^ 1082.9 km/s, a = 0.3361, and ctdmj 
the dark matter velocity dispersion, is related to the 
observed galaxy velocity dispersion by the velocity bias 
parameter, by = ct/ctdm. Simulations generally suggest 
that 0.9 < by < 1.1, and the most recent numerical re- 
sults (Evrard et al. 2008, and references therein) indicate 
that galaxies are essentially unbiased tracers of the dark 
matter potential (i.e., = 1). Adopting this value, we 
obtain M200 = 1-04 ^0 44 x 10^^ Mq. Dynamical masses 
corresponding to the various velocity dispersion measures 
reported in §3.2 are tabulated in Table 2. The Carl- 
berg et al. (1997) relation yields masses a factor of ~ 1.5 
larger. 

Detailed studies of the velocity dispersion profile (e.g., 
Girardi et al. 1993; Fadda et al. 1996) indicate that dis- 



persions measured within radii less than r2oo may be bi- 
ased high. In AlO, raoo for SPT-CL J0546-5345 is esti- 
mated to be « 1.57 Mpc, and we are therefore probing 
an aperture of ~ 0.57r2oo- Numerical (Biviano et al. 
2006) and observational (Rines and Diaferio 2006; Kat- 
gert et al. 1996, 1998) studies suggest that for our early- 
type galaxy sample and the clustcrcentric radius probed, 
the bias is expected to be small. For example, Biviano 
et al. (2006) report the expected bias in the dispersion 
as a function of the fractional virial radius probed; for 
our dispersion measurement the expected upward bias is 
< 2.5%. This implies a corresponding bias in the dy- 
namical mass of < 7%, which would lower the dynamical 
mass to ~ 9.6 x 10^^ Mq. While this may well be a more 
accurate value, the bias correction is both uncertain and 
small. We therefore defer a more extensive analysis to 
a futm-e paper, when a larger sample of spectroscopic 
members out to r2oo is available. 

Finally, we note that no other z > 1 cluster with 
> 10 spcctroscopically confirmed members has a > 1000 
km/s. With a dispersion of a = 1179 km/s, SPT- 
CL J0546-5345 is unambiguously the most dynamically 
massive cluster yet identified at 2; > 1. 

4. DISCUSSION 
4.1. a — Tx Relation 

At Tx = 7.5t\i keV, SPT-CL J0546-5345 is among 

the hottest clusters yet observed in the X-ray at z > 1 
(Rosati et al. 2009). Figure 4 shows the well-known cor- 
relation between X-ray temperature and galaxy velocity 
dispersion (Lubin and Bahcall 1993; Bird et al. 1995; Gi- 
rardi et al. 1996; Horner et al. 1999; Xue and Wu 2000; 
Ortiz-Gil et al. 2004). The clusters shown are the sub- 
set of the well-studied, low redshift X-ray cluster sam- 
ple of Vikhlinin et al. (2009) for which accurate veloc- 
ity dispersion measurements, from at least 40 galaxies, 
are available from Girardi et al. (1996). The temper- 
atures are measured in an identical manner to that of 
SPT-CL J0546-5345, which is plotted as the large red 
circle. The dashed line is the best-fit relation Girardi 
et al. (1996) obtained using previous temperatures from 
the literature, and the solid line is our own fit using 
the Vikhlinin ct al. (2009) temperatures. Both fits ex- 
clude SPT-CL J0546-5345, and are consistent with typi- 
cal measurements (e.g., Lubin and Bahcall 1993; Horner 
et al. 1999). Although the comparison clusters are all 
at low redshift {z < 0.1), there is no expectation of, nor 
evidence for, evolution in this relation to z ~ 1 (e.g., Wu 
et al. 1998; Tran et al. 1999). The measured dispersion 
and X-ray temperature for SPT-CL J0546-5345 faU right 
on the a — Tx relation, offering strong corroboration of 
the very large dynamical mass presented here. 

4.2. X-Ray, SZE and Optical Mass Measures 

We now compare the dynamical mass with several 

other mass measures available for SPT-CL J0546-5345, 
all reported in terms of A'/20o(Pc). A summary of these 
these mass measures is given in Table 2. 

In AlO the cluster mass within rsoo is measured us- 
ing the Yx-M^oo relation (Vikhlinin et al. 2009) to be 
M<ioo,Yx = (5.33 ± 0.62) x lO" Mq. The systematic 
uncertainty of this scaling relation calibration was es- 
timated by Vikhfinin et al. (2009) to be ±9% by com- 
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Fig. 4. — Empirical a — Tx correlation for the X-ray clusters from 
Vikhlinin et al. (2009), with velocity dispersions taken from Girardi 
et al. (1996). The temperatures are measured exactly as for SPT- 
CL J0546-5345. The fits, described in the text, are typical of those 
found by other authors (e.g., Lubin and Bahcall 1993; Horner et al. 
1999). SPT-CL J0546-5345 (red circle) falls right on this relation. 

paring X-ray mass estimates to weak lensing mass mea- 
surements for a representative sample of clusters. We in- 
clude this uncertainty in quadrature when converting to 
M2oo,Fxi s-iid '^^ estimate that M2m,Yx = (8.23 ±1.21) x 
10^^ Mq. This X-ray mass is consistent with the mea- 
sured dynamical mass within the errors. Although the 
X-ray image shows evidence of a possible merger, the 
Yx estimator is expected to be robust to such events 
(Vikhlinin et al. 2009). Indeed, the majority of high 
redshift clusters have asymmetrical X-ray morphologies 
(e.g., Vikhlinin et al. 2009; AID), so the inferred X-ray 
mass of SPT-CL J0546-5345 should be of comparable 
accuracy to other clusters at similar redshifts. We also 
estimate the mass from the temperature, via the Tx- 
M500 relation of Vikhlinin et al. (2009), to be Af2oo,Tx = 
(8.11 ± 1.89) X lOi'^ Mq. 

The SZE measurement of the integrated Comptoniza- 
tion, is expected to be the most accurate SZ mass 
proxy (Motl et al. 2005; Nagai 2006). In AlO, they 
measure a spherically deprojected Comptonization of 
Ysz = (3.5 ± 0.6) X 10" MekeV, where the uncertainty 
is estimated from the quadrature sum of statistical errors 
and the uncertainty in the assumed gas profile. From this 
measurement, we use the Ysz-M^qq relation reported 
in AlO to estimate a mass, where we marginalize over 
the allowed range of normalization and slope parameters 
from their fit. We include an additional ±9% uncertainty 
to account for the systematic uncertainty of the Yx- 
based masses used to calibrate their relation^° Using this 
method, we infer a mass of M2oo,Yg2= (7.19±1.51) x lO^'' 

'"^^ The SZ and X-ray mass estimates aren't completely indepen- 
dent. We expect that their largest correlation will be from the 
mass calibration of the Yx-M^qq and ^52— M50Q relations, which 
are based on the same X-ray measurements from Vikhlinin et al. 
(2009). We add the 9% uncertainty to the Ygz mass estimate to 
account for this, and assume that any other correlated uncertainty 
is negligible. 



A^Q, in agreement with the X-ray and dynamical mass 
estimates. 

Another SZ mass estimate can be inferred from the 
SZE significance-mass scaling relation established via 
WMAP7-constrained simulations in VIO. AlO compared 
this SZE significance-based mass proxy with X-ray in- 
ferred masses for their sample of 15 SZE-selected clus- 
ters, finding that the significance-based SZE masses may 
be biased low, with an average ratio of the SZE to X-ray 
masses of 0.89±0.06. The resuh for SPT-CL J0546-5345, 
M200.SZ = (5.03±1.13±0.77) X IQi'' A/q, where the uncer- 
tainties are statistical followed by systematic, is indeed 
^ 1.6a below the 68% confidence interval of the more 
precise X-ray mass, although it is in agreement with the 
dynamical mass within the errors. Accounting for the 
bias quantified in AlO brings the significance-based mass 
estimate into agreement with the X-ray mass. 

HID present a scaling relation between SZE-derived 
mass and two optical richness measures, iVgai and iV2oo- 
These richness measures count the number of galaxies 
on the red-sequence at the cluster redshift within clus- 
tercentric radii of 1 Mpc and r2ooj respectively. The 
derived masses, A/2oo(Afgai) = (9.2±4.9±2.7) x lO^"* Af© 
and A^20o(A^20o) = (8.5±5.7±2.5) x 10^'' Mq, are consis- 
tent with the independent richness-based mass estimate 
of Menanteau and Hughes (2009) for this cluster. Al- 
though noisy mass proxies, these richness measures offer 
independent evidence of the high mass of SPT-CL J0546- 
5345. 

4.3. Combined Mass for SPT-CL J0546-5345 

Table 2 lists total mass measures for SPT-CL J0546- 
5345 from dynamical. X-ray, SZE and optical richness 
mass proxies. For each physical probe, the bold entries 
indicate the specific proxy expected to yield the most ro- 
bust mass. For dispersions, this is the biweight measure, 
which is optimal for the number of spectroscopic mem- 
bers we have (Beers et al. 1990). For the X-ray, the Yx 
estimator is chosen over the Tx estimator because it is 
expected to be significantly less biased for mergers and 
to have less intrinsic scatter with mass (Kravtsov et al. 
2006; Vikhlinin et al. 2009). For the SZE, the Ysz esti- 
mator is chosen over the S/N estimator because scaling 
relations for it have been calibrated directly from ex- 
ternal observations rather than simulations and it is ex- 
pected to have less intrinsic scatter with mass (Kravtsov 
et al. 2006; VIO). Finally, the A''2oo estimator is the best- 
tested optical richness measure. 

We optimally combine these independent measure- 
ments, weighting them by their errors. We symmetrize 
the errors on the dynamical mass and add the statisti- 
cal and systematic errors in quadrature for the richness. 
The combined mass is insensitive to these details since 
these are the least constrained measures. The resulting 
mass is A^200 = (7.95 ± 0.92) x 10^^ Mq, which we take 
as our best estimate of the mass of SPT-CL J0546-5345. 

4.4. SPT-CL JO546-5345 in a Cosmological Context 

Using the Tinker et al. (2008) mass function we cal- 
culate that SPT-CL J0546-5345, which at z = 1.067 is 
already almost as massive as the Coma cluster, will in- 
crease in mass by a factor of ^ 4 over the next 8 Gyr. 
It should therefore grow into one of the most massive 
clusters in the Universe by the present day. 



SPT-CL J0546-5345: A Massive z>l Galaxy Cluster Selected Via the SZE 
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TABLE 2 

Comparison of Mass Measurements for SPT-CL J0546-5345 











Mass Scaling 




Mass Type 


Proxy 


Measurement 


Units 


Relation 


(10" Mq) 


Dispersion 


Biweight 


1 1 7Q +2'*2 


km/s 


cr- Af200 (Evrard ct al. 2008) 






Gapper 


11701?^° 


km/s 


0--M200 (Evrard et al. 2008) 


10.1 las 




Std Deviation 


11381?°^ 


km/s 


C-M2Q0 (Evrard et al. 2008) 


q o+B.b 


X-ray 


Yx 


5.3 ± 1.0 


xlO" M0keV 


Yx-Mr,QQ (Vikhlinin et al. 2009) 


8.23 ± 1.21 




Tx 


7.511:1 


keV 


TX-M500 (Vikhlinin et al. 2009) 


8.11 ± 1.89 


SZE 




3.5 ±0.6 


xlO" MgkeV 


^SZ - ^500 (AlO) 


7.19 ± 1.51 




S/N at 150 GHz 


7.69 




i - M500 (VIO) 


5.03 ±1.13 ±0.77 


Richness 


^200 


80 ±31 


galaxies 


JV200 - M2Q0 (HIO) 


8.5 ±5.7 ±2.5 






66 ±7 


galaxies 


ATgal - M200 (HIO) 


9.2 ±4.9 ±2.7 



Best 



Combined 



7.95 ± 0.92 



^ M500 masses were scaled to M200 masses assuming an NFW density profile and the mass-concentration 
relation of Duffy et al. (2008). 

^ We do not correct the dynamical masses for the small potential bias in the dispersion discussed in 
§3.4. 



Historically, the existence of individual massive galaxy 
clusters at high redshift has been used to constrain cos- 
mological models (e.g., Donahue et al. 1998). While we 
are entering an era where large samples of SZ clusters will 
soon be available in addressing such questions, it is in- 
teresting to ask whether the existence of SPT-CL J0546- 
5345 is expected in our survey volume for a concordance 
ACDM cosmology. Convolving the Tinker ct al. (2008) 
mass function with a Gaussian mass probability func- 
tion for the best-estimate mass of SPT-CL J0546-5345, 
we find that ACDM predicts the existence of 0.18 clusters 
of this mass or higher at 2: > 1 in 178 deg^. We note that 
the M500 — Yx relation was calibrated at lower redshifts 
(z < 0.6) than SPT-CL J0546-5345, and the accuracy of 
this relation has not been verified at z > 1. Given the 
caveats and uncertainties inherent in this calculation, we 
conclude that the existence of SPT-CL J0546-5345 in our 
survey volume is unsurprising. The complete SPT sam- 
ple, combined with improved mass scaling relations based 
on complementary mass measures, will permit ACDM 
predictions for the high end of the cluster mass function 
to be robustly tested at high redshift. 

5. CONCLUSIONS 

We report the spectroscopic confirmation SPT- 
CL J0546-5345 &t z = 1.067, the first SZE-selected 
galaxy cluster at 2; > 1. We measure a robust ve- 
locity dispersion from 18 early-type members of cr = 
1179 l^g7 km/s, corresponding to a dynamical mass of 



M200 = 1.0 



-1-0.6 
-0.4 



X 1015 Mq. SPT-CL J0546-5345 is the 
most dynamically massive cluster yet identified, from any 
method, at z > 1. 

We find excellent consistency across several indepen- 
dent mass measures. The measured dispersion and X-ray 
temperature of SPT-CL J0546-5345 fall right on the a- 
Tx relation. The X-ray, SZE and richness-based mass 
estimates are all consistent with the dynamical mass. 



and with each other, within the errors. Combining all 
the mass measures, we derive a best-estimate mass for 
SPT-CL J0546-5345 of M200 = (7.95 ±0.92) x lO" Mq. 

In ACDM we expect 0.18 clusters consistent with this 
mass above z > 1 in our survey area. Given the un- 
certainties in X-ray scaling relations at high redshift, we 
conclude that the existence of SPT-CL J0546-5345 in 
our survey volume is unsurprising. The complete SPT 
sample will provide the large sample of clusters required 
to robustly test the high redshift, high mass end of the 
cluster mass function. 
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